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Abstract. We describe a simple technique for optimizing the extraction of the forward-backward asym- 
metry {Aft) of Drell-Yan lepton pairs (e + e~, n + fi~) produced in pp and pp collisions at hadron colliders. 
The method employs simple event weights which are functions of the rapidity and | cos#| decay angle of 
the lepton pair. It yields the best estimate of the acceptance corrected parton level (qq) forward backward 
asymmetry as a function of final state dilepton mass (Mu). Typically, when compared to the simple count 
method, the technique reduces the statistical errors by 20% for pp, and 40% for pp collisions, respectively. 
The method is equivalent to using the maximum likelihood method, but is much easier to implement. The 
technique can be used to search for new high mass and large width Z' bosons which may be best detected 
through the observation of deviations from the Standard Model expectation for the forward-backward 
asymmetry. The technique can also be applied in the extraction of the toward- barckwad asymmetry in the 
production of top-antitop pairs. 
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1 Introduction 

The Drell-Yan process in which qq annihilations form in- 
termediate 7* or Z (7* /Z) vector bosons decaying to lep- 
ton (e + e~, pairs is particularly useful in search- 
ing for new interactions at large momentum transfers, 
Q 2 = M 2 e , where Mu is the invariant mass of the lepton 
pair. In leading order (LO) approximation, the momen- 
tum fractions X\, X2 carried by the initial state quarks 
and antiquarks in the proton and antiproton/proton, re- 
spectively, are related to the rapidity \y\ pQ of the 7* jZ 
boson via the equation x\^ — (Mu/ ' 1 /s)e ± ' y ' , where ^/s is 
the center of mass energy. Dilepton pairs produced at large 
\y I originate from collisions in which one parton carries a 
large and the other a small momentum fraction x. 

Drell-Yan lepton pairs which are produced in qq anni- 
hilations display a forward-backward asymmetry because 
of the interference between photon and Z boson exchange[2] 
This forward-backward asymmetry would be modified by 
new resonances (e.g. additional heavier Z' bosons [2J) or 
new interactions at large mass scales. 

Although the mass limits from LEP [3] and Tevatron 
(pp) experiments CDF [4] and and D0 [5] for a variety of 
Z' models are in the 0.6-0.8 TeV/c 2 range, the limits are 
much lower if the Z' width (typically Ez' ~ 0.01 • Mz>) is 
increased to account for the possibility of additional de- 



cays modes to exotic fermions (which are predicted in Eg 
models [6]), and/or supersymmetric particles. The limits 
are even lower if one includes the possibility of a more 
general model with enhanced couplings to the third gen- 
eration. 

Such a Z' (which has larger width e.g. Tz> =0.1 Mz<) 
would produce only a small signal in the dilepton mass 
spectrum because the total cross section is proportional 
to the square of the amplitude. However, the change in 
the forward-backward asymmetry which results from the 
interference with the standard model process is linearly 
proportional to the amplitude and would be observable 
as a change in the forward backward asymmetry. This 
change will occur around the mass of the Z' boson, and 
also in some mass range below and above the mass of the 
Z' boson [7J. This point is illustrated in Figure 1 which is 
taken from Reference jTJ. The figure compares early (120 
pb~ x ) high mass Drell-Yan da/dMu data from CDF and 
DO, and early Afb data from CDF to Standard Model 
theoretical predictions and (as an example) to a prediction 
with an extra Eq boson with Mz< = 350 GeV/c 2 and 
T Z ' = O.lMz', for <j> = 60° (solid) and = 173° (dotted). 
As can be seen for this case, the signal in Afb is larger than 
the signal in da/dMu- Since such new particles or new 
interactions may be best detected through the observation 
of deviations from the Standard Model expectation for 
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small statistical samples. The method is equivalent to us- 
ing the maximum likelihood method, but is much easier 
to implement. The technique can also be applied in the 
extraction of the foward-backward asymmetry in the pro- 
duction of top-antitop pairs. 



2 qq annihilations 

The differential cross-section for the parton level process 
for qq annihilation can be written as 
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where 9 is the emission angle of the positive lepton rela- 
tive to the quark momentum in the center of mass frame, 
and A and B' are parameters that depend on the weak 
isospin and charge of the incoming fermions. (We use B' 
for coefficient of the cosO at the parton level, and B for 
the coefficient of the cos9 at the nucleon level). 

The cross sections for forward (for—) events (ap) and 
backward (back—) events (&b) are given by 



1. Early Tevatron data (120 pb' 1 ). (a) da/dM u for 
e + e~ pairs from CDF and DO, and \i~ pairs from CDF. 
The Standard Model theoretical predictions (dashed line) have 
been normalized (by a factor of 1.11) to the CDF data in the 
Z boson mass region, (b) CDF data for Afb versus Ma (e + e~ 
pairs). Also shown is the standard model expectation (dashed 
line) and the predictions for dcr/dM and Afb with an extra E§ 
boson with M z , = 350 GeV/c 2 and r z , = 0AM z >, for 4> = 60° 
(solid) and <f> = 173° (dotted). 



the forward-backward asymmetry, it is useful to devise 
experimental techniques to measure the forward-backward 
asymmetry (Afb) with the best possible precision. Since 
fine mass bins (~ 25GeV/c 2 ) are required, the number of 
events per bin at large Mm is small and the measurements 
are statistically limited. 

In this communication we describe a simple technique 
for optimizing the extraction of the forward-backward asym- 
metry in the production of Drell-Yan dilepton paris. The 
method employs simple event weights, which are functions 
of the rapidity and | cos6>| decay angle of the lepton pair. 
The method yields the smallest statistical uncertainty in 
the measurement of the forward backward asymmetry as 
a function of Ma- It can be directly applied to current pp 
data at the Fermilab Tevatron, as well as new data that 
will be collected in pp collisions at the Large Hadron Col- 
lider (LHC). The method works well for both large and 
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The electroweak interaction introduces the asymmetry 
(a linear dependence on cos 9), which can be expressed as 
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For pp collisions (e.g. at the Tevatron), the direction 
of the quark is predominately in the proton direction, and 
the direction of the antiquark is predominately in the an- 
tiproton direction. Therefore, the forward backward asym- 
metry for qq processes is easy to measure in pp collisions. 

If Nf is in number of events in the forward direction of 
the quark and is the number of events in the backward 
direction of the quark we obtain the following expression 
for the total forward backward-asymmetry (Aj° b tal ) and its 

error (AAf b tal ): 
[AfA 
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QCD corrections also dilute the forward-backward asym- 
metry as follows: 
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where we have used AN f = (iV/) 1 / 2 and AN b = (N b ] 
and N — Nf+N b . Since for Poisson statistics [TT], the frac- 
tional error is (l/7V e a:pected) 1/2 and not (l/AU se ™ et 2) 1/2 , 
we use A f b(expected) m equation 7j For pp collisions above 



the Z mass peak, A 



fb(expected) " 



j. In this region, AA 



fb 



0.800 ■ (1/A0 1 / 2 . 

Therefore, a measurement with 100 events yields a sta- 
tistical error of 0.08. This level of precision is needed to 
observe the deviation from the Standard Model for the Z' 
example shown in figure 1. Later in this paper we show 
that a reduction in the error (of about 20%) can be ob- 
tained by using the information in the angular distribution 
of the forward and backward events. 
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On average, both QCD correction terms are very small. 
For Drell-Yan events the average ratio Pt/Mu is about 
(7) 0.13 (e.g. 12 GeV P T at the Z mass of 90 GeV). This im- 
plies that on average < A (P T ) >= 0.018, and < jfr(Pr >= 
0.991. 

V 2 The emission of photons by the final state leptons also 
distorts the angular distribution (in addition to smear- 
ing events to lower mass) . The effects of final state QED 
photons are best accounted for as part of an overall radia- 
tive correction. Alternatively, the angular distortion from 
QED radiation can also be corrected for by using effective 
angular coefficients which combine the QCD and QED 
corrections (expressions for the larger A^ ad (PT, M«) and 

larger ^-(Py, Mil) are given in Appendix A). 



4 Misidentification of the quark direction in 
pp collisions 



3 Quarks bound in a nucleon 

When quarks are bound in the nucleon, the dilepton can 
be produced with non-zero transverse momentum. This is 
described in detail in appendix A. For pp or pp collisions 
the angular distribution of 7*/Z vector bosons decaying 
to e + e~ or fi + p~ pairs is given by: 

- d ° - = All + cos 2 6 + Mff)] +Bcos6 (8) 
d(cosu) 

h(6) = ^A (M u ,P T )(l~3cos 2 9) (9) 

The qq center of mass frame is well defined when the lep- 
ton pair has zero transverse momentum (Pt)- For a non- 
zero transverse momentum of the dilepton pair, the qq 
center of mass frame is approximated by the Collins-Soper 
frame [8]. 

The term h(6, Ma, Pt) is a small QCD correction term 
which is zero when the transverse momentum of the dilep- 
ton pair is zero. As shown in Appendix A, the angular coef- 
ficient Aq is only a function of the dilepton mass (Mu) and 
transverse momentum (Pt). The h(0, Ma, Pt) term inte- 
grates to zero when the cross section is integrated over all 
cos 2 9. and increases with transverse momentum. For sim- 
plicity we will use A = (e.g. Pr=0) in our calculations. 
In order to include the dependence on the h(6, Mu, Pt) 
term, we can use the following function derived in Ap- 
pendix A: 

A ° 1 + P 2 /M 2 (10) 



For pp collisions the direction of the quark is primarily 
along the direction of the proton. However, there is a 
a small probability for a misidentification (misID) that 
originates from the charge misidentification probability in 
the tracker. There is also an additional small misID that 
comes from the small fraction of events in which a sea 
antiquark in the proton interacts with a sea quark in the 
antiproton. This misID is very small for large dilepton 
final state mass. 

We now turn to the measurement of the forward-backward 
asymmetry in pp collisions (e.g. at the Large Hadron Col- 
lider). Although pp collisions are symmetric, there is still 
a forward backward asymmetry if the quark direction is 
defined to be the direction of motion of the Drell-Yan pair. 
This originates from the fact that on average, quarks carry 
a larger fraction of the proton's momentum than anti- 
quarks. However, there is a significant misID that origi- 
nates from the fraction of events for which the antiquark 
carries a larger fraction of the momentum than the quark. 
This misID dilutes the observed asymmetry. 

In most theoretical studies of the production of new Z' 
bosons inpp collisions at the LHC, this dilution is included 
in calculation for the prediction for the observed forward- 
backward asymmetry. Here, we show that we can obtain 
a higher sensitivity to new particle searches by correcting 
the data for the misID fraction on an event by event 
basis. 

At small rapidity the misidentification probability Wi 
is large (0.5 at \y\ = 0). At large rapidity the misidenti- 
fication probability Wi is small. We show below that by 
taking this information into account we can reduce the 
error on the extracted qq asymmetry. 
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5 Correcting for misID in pp collisions 

We illustrate this point for the case of a high statistics 
measurement of the forward backward asymmetry in a 
specific mass bin (e.g 450 GeV/c 2 ). 

We first extract the qq asymmetry by correcting the 
measured asymmetry for the average misID probability 
(this commonly used method is called the event count 
method). 

We then show that we reduce the error on the ex- 
tracted quark-antiquark asymmetry by binning the data 
in ten rapidity bins and fitting for the weighted average 
of the extracted parton level asymmetries from all of the 
ten rapidity bins. 

Then we derive an event weighting technique that is 
equivalent to the fit method in the high statistics case, but 
which can also be used in the limit of very small statistical 
samples. 

For proton-antiproton collisions we define the misID 
probability (wi) for each yi bin as the fraction of events 
for which the antiquark carries a larger fraction of the 
momentum than the quark. (For proton-antiproton colli- 
sions Wi is the fraction of events for which the interactions 
occurs between a sea quark in the antiproton and a sea 
antiquark in the proton). 



M e , 



51(2) 



x e 



±lwl 



w, 



avor 

v q {q(x 2 ) ■ q(xi)} 



E flavor VqUfa) ' + ?Ol) ' q(x 2 )} 



(12) 



(13) 



Where M« is the dilepton mass (M^ or M ee ), and 
is the rapidity of the dilepton pair. Here q(x) denotes the 
quark distributions (u(x), d(x) , s(x) , c(x),_b(x)) and q(x) 
denotes the antiquark distributions (u(x), d(x), s(x), c(x), 
b(x)) for the various flavors in the nucleon. The parameter 
v q denotes the Z/7 couplings of to each flavor (which a 
function of the dilepton mass). At large M« the u(x) and 
u(x) quark distributions dominate the expression for the 
dilution factor. 

Note that the misID in the forward-backward asym- 
metry should include the interference between photon and 
Z boson exchange for each quark flavor, which in general is 
a more complicated function of the couplings. In addition, 
the misID is affected by radiative emission of photons and 
detector resolution. Therefore, it is best to use a Monte 
Carlo generator (such as PYTHIA[9] or ZGRAD2[10]) to 
empirically determine the functional dependence of the 
dilution factor w, as a function of the two measured vari- 
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We now proceed to correct for the misID and extract 
the qq forward-backward asymmetry for each one of the 
ten yi rapidity bins. In the expressions below, rif^ and nb,i 
are the defined as the measured (i.e. diluted) number of 
forward events and backward events in each bin, and Nf i 
and Nb i are defined as the number of true forward and 
true backward events (for qq collisions) in the bin. 



For a given yi rapidity bin with a misidentification 
probability W4 the measured and true number of forward 
and backward events are related by the following expres- 



nf,i = Nf >i (l-Wi) + N b>i (w i ) (15) 
n b j = N bji (l - Wi) + N {ti (wi) 
N f,i = n fA l - Wi)/Li - n b ,i(wi)/Li 
Nbj = - Wi)/Li - nf ti (wi)/Li 

where Li = (1 — 2w{) is defined as the dilution factor. 
For a \yi\ rapidity bin, the corrected parton level asym- 
metry is given by 
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where we have used Arif i = (rif i) 1 / 2 and An o i — (n^ i) 1 / 2 . 

We find that a measurement of the qq asymmetry in 
the case where there is misID probability of Wi results 
in an increase of the error in the extracted parton level 
asymmetry by a factor of I /Li (which is equivalent to 
reducing the number of events by a factor of L 2 = (1 — 
2 Wl ) 2 ). 

If we want to combine different \yi\ bins together, we 
need to weight the events by the inverse of the square of 
the statistical error in each bin. This is achieved by multi- 
plying the expressions for Nfj and N bii by L 2 . Since this 
factor appears both in the numerator and denominator of 
the expression for Aj ^ , it does not change the extracted 
value or error of the parton level asymmetry. However, 
when we combine \yi\ bins together using event weighting, 
this factor accounts for the difference in statistical errors 
between the \yi\ bins as follows. 



ki.i = (1 - Wi)(l - 2wi) 
&2,i = (Wi)(l - 2Wi) 

Ntotal = J2 W 
all— events 

S f= H Ki - H fc 2,i 
for— events back— events 

[AS f f= £ kl+ £ kl 

for— events back— events 

Sb = ^1,2 - ^ 2 '* 

back— events for— events 



(19) 
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\As b ] 2 = E k h + E k h 

back — events f or— events 

a total _ g/ ~ ^fr 

A f» ~ S f + Sb 

Now [/iS/] and [ASb] are correlated with each other in 
a complicated way. In order to simplify the calculation of 
the error, we combine terms to isolate sums which are for 
forward events, and sums which are for backward events, 
separately as follows: 
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Now AAi is 100% correlated with AB\ and AA 2 is 
100% correlated with AB 2 . We handle these correlations 
as follows. 
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A specific numerical example is shown in Table [T] Here 
we show the case of a measurement of Afb in ten bins of 



\y\ for the range of \y\ between and 2.0. For this study 
we assume that an asymmetry of Afb = 0.6 is measured 
with 1000 events (for < |y| < 2) in pp collisions at 
the LHC. We assume that the differential cross section is 
constant in \y\. The misID probability is assumed to be 
ui(|y|) = 0.5(2— |y|)/2 which is a simple approximation to 
to the misID values for a dilepton mass of 450 GeV/c 2 
at the LHC. In each range in \y\ we compare the error 
in the extracted pp asymmetry from a simple count of 
events; the error extracted from a least square fit to the 
values extracted from each \yi\ bin, and the error from 
the our event weighting formula. As expected, the error 
from the least square fit to the ten bins is the same as 
the error from our weighting formula. The error using a 
simple count is about 20% larger than the error using the 
event weighting scheme. The last column shows the error 
for the case oi pp collisions (with a misID=0). 

Our method is equivalent to combining different \yi\ 
bins together by weighting the events by the inverse of the 
square of the " expected" statistical error in each bin. This 
has been done by multiplying the expressions for Nf i and 
Nb,i by Lf. Weighting by a statistical error which is based 
on the "expected" number of events in a bin (instead of the 
measured number of events in a bin) is more correct, and is 
equivalent [TT] to using the maximum likelihood technique 
for combining bins. Therefore, our event weighting tech- 
nique is equivalent to using a maximum likelihood tech- 
nique, but is much easier to implement. 



6 Including information in the angular 
distribution in pp collisions 



We now investigate how much can be gained by looking at 
the asymmetry in bins of \c 3 ■] — \ cos 8j \ . We start with the 
case of qq collisions and divide the sample into ten bins in 
| cos9j\. The asymmetry as a function of \cj \ bin is: 



°f(\ c j)\ +°b\{Cj\) 
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At \cj\ = | cos ^ | = 0, the measured asymmetry A/b_j(0)=0. 
| cos 0j | = 0.45, the measured asymmetry A fb-j (0.45) = 
= |cos#j| = 1, the asymmetry Afb(l) — 



At , , 

Atotal Ai 



(4/3)^4^*°'. The measured asymmetry in each \cj \ — cos 9j 
bin can be related to the total (integrated over all cos6>) 
asymmetry. Therefore each | cos 9j \ bin provides an inde- 
pendent measurement Ajf^j of the total asymmetry. 
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Table 1. Proton-Proton collisions: Numerical example of measuring Aft, in ten bins of \y\ for the range of \y\ between 
and 2.0. For this study we assume an asymmetry Af b — 0.6 which is measured with 1000 events (for < \y\ < 2) in pp 
collisions at the LHC. We assume that the differential cross section is constant in \y\. The misID probability is assumed to be 
w(|j/|) = 0.5(2 — |j/|)/2 which is a simple approximation to to the misID values for a dilepton mass of 450 GeV/c 2 . In each 
range in \y\ we compare the error in the extracted qq asymmetry from a simple count of events, the error using a least square 
fit to the values extracted for all bins, and the error from the proposed event weighting formula. The error using a simple 
count is about 20% larger than the error using the the proposed event weighting scheme. The last column shows the error for 
the case of pp collisions (with a misID=0). 



\y\- range 


N 


simple Count 


\y\ bins fit 


Weights 


Improv. 


pp 


range 


bins 


Error pp 


Error pp 


Error pp 


Factor pp 


Error 


0.-0.2 


1 


1.999 


1.9991 


1.9991 


1.00 


0.087 


0-0.4 


2 


0.706 


0.6301 


0.6298 


1.12 


0.061 


0-0.6 


3 


0.383 


0.3350 


0.3346 


1.14 


0.050 


0-0.8 


4 


0.248 


0.2145 


0.2141 


1.15 


0.043 


0-1.0 


5 


0.177 


0.1512 


0.1511 


1.16 


0.039 


0-1.2 


6 


0.134 


0.1136 


0.1131 


1.17 


0.035 


0-1.4 


7 


0.106 


0.0886 


0.0882 


1.18 


0.033 


0-1.6 


8 


0.086 


0.0711 


0.0706 


1.21 


0.031 


0-1.8 


9 


0.071 


0.0582 


0.0578 


1.23 


0.029 


0-2.0 


10 


0.060 


0.0483 


0.0479 


1.26 


0.027 



Table 2. Proton- Antiproton collisions: Numerical example for an asymmetry Aft — 0.6 measured with 1000 events. Here Af° b tal 
is measured in ten bins of | cos#|. In each | cos8j\ bin we compare the error from the standard error formula, and the error from 
the event weighting formula. As expected, the two yield identical results. In addition, we show a comparison of the average 
j^otai Q £ a jj ^ en | cog | kjjjg calculated two different ways. The error in the average extracted from a least square fit to the 
10 Af^ values (0.0196) is close to the error in the average determined from the weighted sum of all the events (0.0210). The 
error in Aj° h ial from the weighted sum of all the events is 20% lower than the error of 0.0253 obtained from a simple count of all 
forward and backward events. 



\Ci\ = 

cos 9j 


Uf 


n b 


Afb-j{\ci\) 
measured 


A total 
A fb-j 

extracted 


A a total 

| cos 9j | bin 


A A total 
nA fb-j 

from weights 


0.05 
0.15 
0.25 
0.35 
0.45 
0.55 
0.65 
0.75 
0.85 
0.95 


41 
47 
55 
63 
72 
82 
92 
104 
116 
128 


34 
29 
25 
21 
18 
16 
14 
14 
14 
14 


0.080 
0.235 
0.376 
0.499 
0.599 
0.676 
0.731 
0.768 
0.790 
0.799 


0.60 
0.60 
0.60 
0.60 
0.60 
0.60 
0.60 
0.60 
0.60 
0.60 


0.864 
0.284 
0.165 
0.114 
0.085 
0.066 
0.054 
0.045 
0.041 
0.038 


0.864 
0.284 
0.165 
0.114 
0.085 
0.066 
0.054 
0.046 
0.041 
0.028 


all 


800 


200 


0.6 


0.60 


0.0196 


0.0210 


all 


800 


200 


0.6 


0.60 


simple-count 


0.0253 



a Atotal 
nA fb-j 



8Mj N fij + Nt 



N f,j N b,j 



1/2 



M, = 



(l + c? + h(M u ,6,P T )) 

where we have used ANf j = (Nf j ) 1 / 2 and AN b j = 

The above expression shows that for case in which we 
have same number of events in each of the | cos 0j | bins, 
the error in the extracted measurement of A*^_j from the 



data in a specific 



| cos 6j | bin is equal to the error 



of the measured asymmetry in the bin divided by a factor 



— (i+c?+h(6)) ' f actor comes from the fact that 

the extracted total forward-backward asymmetry is more 
sensitive to events at large \cj\ = | cos0j|. 

We now convert the procedure to event weight tech- 
nique. We define Af^ = {3/8)(N A , J /N BJ ). 



N f , = N Ad (l + c 2 + h(6)) + \ Cj \N BJ (24) 



N h ,=N A ,{l + c' + h(8))-\c 3 \N B , 
From which we get: 



+ 



2(1 + c 2 + h) 2(1 + c 2 + h) 



(25) 
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Table 3. Proton- Antiproton collisions: Numerical example of measuring Afb for different acceptance ranges in |cos#|. The 
example is for an asymmetry Afb — 0.6 measured with a total of 1000 events (for all values of | cos#|). For each | cos#| range 
we show the error in Afb from a simple count, the error from the event weighting procedure, and the error from fitting bins in 
| cos #|| (which is very similar to event weighting). Also shown is the improvement factor in the error when the event weighting 
procedure is used (versus a simple count). For a typical range of | cos#|, using the simple weighting formula leads to more than 
20% reduction in the error. 



| cos 6»| 


N 


simple Count 


| cos 6\ bins 


Event Weights 


Improvement 


range 


bins 


Error 


fit Error 


Error 


factor 


0-0.1 


1 


0.8642 


0.8642 


0.8642 


1.000 


0-0.2 


2 


0.3042 


0.2796 


0.2687 


1.132 


0-0.3 


3 


0.1644 


0.1410 


0.1403 


1.172 


0-0.4 


4 


0.1058 


0.0884 


0.0881 


1.201 


0-0.5 


5 


0.0749 


0.0611 


0.0611 


1.225 


0-0.6 


6 


0.0563 


0.0459 


0.0454 


1.241 


0-0.7 


7 


0.0442 


0.0346 


0.0354 


1.247 


0-0.8 


8 


0.0358 


0.0277 


0.0288 


1.242 


0-0.9 


9 


0.0298 


0.0230 


0.0242 


1.227 


0-1.0 


10 


0.0253 


0.0196 


0.0210 


1.205 



N B ,j = 



N 



./■:; 



2 k 



2k 



In order to properly weight events for different | cos 9j \ 
bins by the inverse of the square of the error for each bin 
we multiply the above expressions by 



and get : 



Mf 



(1 



cj + h(M u ,e,p T )y 



z i,j = 



^,3 



1 

2 (1 + c§ + /i(0)) 3 

1 kl 

2 (1 + c] + h(9)Y 



A 2 = N b ,j ■ (z ld ) 
B x = N fd ■ (z 2 , j ) 
B 2 = N bd ■ (z 2lj ) 



[AAx] 2 = N 



[AA 2 



N, 
N 



b,j ' z l,j 
fJ ' Z 2,j 



[AB 2 ] 

Aj = Nfj(z hj ) + N bd (z hj ) 
= Ai + A 2 



B, 



Nf,j(z2,j) -N bij (z2 d ) 



= Bi — B 2 



A total 
A fb-j 



3B, 



ZB 1 -B 2 



Aj 8 Ai + A 2 



(26) 



(27) 



(28) 
(29) 



Now AA X is 100% correlated with AB X and AA 2 is 
100% correlated with AB 2 . We handle these correlations 



as follows. 



AA X = ABi 



AA 9 = AB 9 ■ 



Al 

m 

A2 

m 



(30) 



1 



(Ai+A 2 ) 



[El + R 



Eo = 



l^±f(A2Bl + AlB2) 2 
B\ 



[AB 2 

R2 
D 1 



-{A2B\ + A1B2) 2 



Table[2]shows the results of a numerical example for an 
asymmetry Af b = 0.6 measured with 1000 events (we as- 
sume h=0). Here A^° b tal is measured in ten bins of | cos0|. 
In each | cos 0j \ bin we compare the error from the stan- 
dard error formula, and the error from the event weighting 
formula. As expected, the two methods yield identical re- 
sults. 

In addition, we show a comparison of the average Ay^ 0,1 
for the ten | cos 9j \ bins calculated in three different ways. 
The error in the average extracted from a least square 
fit to the 10 Af^ values (0.0196) is close to the error 
in the average determined from the weighted sum of all 
the events (0.0210). The error in 

A total from thc 

weighted 

sum of all the events is 20% lower than the error of 0.0253 
obtained from a simple count of all forward and backward 
events. 

For the case of low statistics, we can use the event 
weighting technique to combine all the events at all value 
of | cos 9 1 (we do not need to bin the events in |cos#|). 
The following are the formulae to extract the best value 
and error from the entire range in | cos#| using the event 
weighting technique : 
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z l,j 
z 2,j 

N total 

At 

A 2 

B x 

B 2 
[AAA 2 = 
[AA 2 ] 2 = 
[ABA 2 = 



[AB 2 



[Af, 



A 
B 

-\total 



[AAf b tal ] 



E 2 



E 



2(1 + £ + h{6,P T )Y 



2{l + c?. + h{6,P T )f 

E w 



all — evejits 

E 

forward— events 



ick — events 



E M 

f or ward— events 

E fa] 

back — events 

E [*U 

f or ward— events 

l~~2 1 



E 

back — events 



E 

f orward— events 

E VU 

back — events 

Ai +A 2 
Bi — B 2 

SB 3 B\ — B 2 
8l ~ 8 M + A 2 



AAi = AB X 
AA 2 



ABo 



Al 

m 

A2 
B2 



(Ax+AzY 



[El + El] 



\- AB £ {A2BI + AIB2) 7 



B\ 
[AB 2 

R2 

u 2 



-(A2BI + A1B2) 2 



(31) 



Note that since we add up the forward and backwards 
events in separate sums, the weighting factors Zij and z 2 j 
are functions of the absolute value | cos#|. 

Table [3] shows a numerical example of measuring Afb 
for different acceptance ranges in | cos 9\ (we assume h=0). 
The example is for an asymmetry Afb = 0.6 and an a sam- 
ple of 1000 events (for < |cos#| < 1). For each range 
of acceptance in | cos#| we show the error in Afb from a 
simple count, the error from the event weighting proce- 
dure, and the error from fitting bins in |cos#| (which is 
very similar to event weighting). Also shown is the im- 
provement factor in the error when the event weighting 
procedure is used (versus a simple count). For a typical 



range of acceptance in | cos#|, using the simple weighting 
formula leads to more than 20% reduction in the error. 

Note that when we use the angular distribution weights, 
the extracted A^° h tal is automatically corrected for the ac- 
ceptance in | cos 9 1 (since the acceptance cancels to first 

order). However, N to tai = Y^aii-events M is ecmal to tne 
observed number of events and is not corrected for accep- 
tance. For Ntotai j an acceptance correction is needed for 
the determination of da/dMu 



7 Combining misID weighting and angular 
distribution weighting in pp collisions 

In pp collisions each event is can be characterized by a 
misID factor which is related to the quark and 

antiquark distribution {x\^ = (M u /y r s)e ± ^) at its value 
of \yi\. In addition, each event has a measured value of 
\cj\ = | cos9j\ . The expressions for combining events with 
different \cj\ = \ cos9j \ and misID Wi values are given by: 



k A,i 
Ntotai 

= 



fci.i - k 2 ,i = (1 - 2wi) 2 
h,i + k i,% = (1 - 2wi) 

E w 



(32) 



all — events 



Z 2,J 

Ai 
A 2 
Bi 
B 2 
[AAA 2 
[AA 2 ] 2 
[ABA 2 

[AB 2 ] 2 

A 
B 

[AAA 
[AA 2 ] 



2 (1 + <?j + h(6)) 3 

i M 

2(l + c 2 j +h{6)) 2 
E [ z hi k A,j] 

f or— events 

E i z hi k A,j} 

back— events 

E i z 2,jk B ,j] 

for— events 

E [ Z 2j fc Bj] 
back— events 

E [ z tj k %j\ 

for— events 

E 

back— events 

E i z h k B,j] 

for— events 

E [ Z 2,j k B,j\ 
back— events 

Ai + A 2 
B\ — B 2 

r „ , Al 

[ABA 



[AB 2 



BI 
A2 

B2 
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Table 4. Proton-Proton Collisions: Numerical example of measuring Afi integrating over various ranges of | cos#| and various 
ranges of \y\ for the case of a dilepton mass of 450 GeV/c 2 at the LHC. Here we assume that the asymmetry Aft = 0.6 is 
measured with 10 6 events for values of |cos#| and < \y\ < 2. Shown is the improvement factor in the the error (over the 
simple count method) when we use event weighting in both | cos#| and w — MisID as a function of \y\. For a typical range of 
acceptance in j cos#| and \y\, the weighting formula leads to a 40% reduction in the error (over the error obtained from a simple 
count). 



y — range 


0-0.2 


0-0.4 


0-0.6 


0-0.8 


0-1.0 


0-1.2 


0-1.4 


0-1.6 


0-1.8 


0-2.0 


j cos 6<| 






















range 






















0-0.1 


1.00 


1.12 


1.14 


1.15 


1.16 


1.17 


1.18 


1.21 


1.23 


1.26 


0-0.2 


1.13 


1.27 


1.29 


1.30 


1.31 


1.32 


1.34 


1.37 


1.39 


1.43 


0-0.3 


1.17 


1.31 


1.34 


1.35 


1.36 


1.37 


1.38 


1.42 


1.44 


1.48 


0-0.4 


1.20 


1.35 


1.37 


1.38 


1.39 


1.41 


1.42 


1.45 


1.48 


1.51 


0-0.5 


1.23 


1.37 


1.40 


1.41 


1.42 


1.43 


1.45 


1.48 


1.51 


1.54 


0-0.6 


1.24 


1.39 


1.41 


1.43 


1.44 


1.45 


1.46 


1.50 


1.53 


1.56 


0-0.7 


1.25 


1.40 


1.42 


1.43 


1.45 


1.46 


1.47 


1.51 


1.53 


1.57 


0-0.8 


1.24 


1.39 


1.42 


1.43 


1.44 


1.45 


1.47 


1.50 


1.53 


1.56 


0-0.9 


1.23 


1.37 


1.40 


1.41 


1.42 


1.44 


1.45 


1.49 


1.51 


1.55 


0-1.0 


1.21 


1.35 


1.37 


1.39 


1.40 


1.41 


1.42 


1.46 


1.48 


1.52 



a total 
A fb 



[AAf tal 



3B 
8 A 
3 



3 B l - B 2 



%A 1 
1 



A? 



{Ax+A 2 f 



[El + El] 



El 



[Agi 
Bl 

[AB 2 . 
ft 2 



\A1B\ + A1B2) 2 



(A2B1 + A\B2f 



Table 4] shows a numerical example of the improve- 
ment in the errors that that can be realized by using the 
information for both the mis ID (wi) and | cos ^ | on an 
event by event basis, versus the measurement of Afb using 
a simple count of all events for all values of | cos 9\ and \y\ 
(within the experimental acceptance). We show the case 
for an asymmetry Afb = 0.6 measured with 10 6 pp events 
(over all values of |cos#| and < \y\ < 2 ) for the case 
of a dilepton mass of 450 GeV/c 2 at the LHC. The two 
dimensional table shows the improvement factor in the 
the errors (over the simple count method) when we use 
event weighting in both | cos#| and MisID as a function 
of \y\. For a typical experimental acceptance in | cos#| and 
|y|, the weighting formula leads to a 40% reduction in the 
error (over the error obtained from a simple count). 

The following are the advantages for using the event 
weighting. 

1 . For pp or pp collisions, if event weighting is done only in 
| cos 9\, the error in the asymmetry is typically reduced 
by a factor of 1.2. This the method that should be 
used to extract the quark couplings from the average 
integrated asymmetry over all values of \y\. 

2. For pp collisions, if the event weighting is be done in 
both | cos£?| and \y\, the error in the asymmetry is re- 
duced by a factor of 1.4. For pp collisions if the weight- 
ing is done in both | cos 8\ and \y\, there is also a further 
reduction in the error. 



3. For pp or pp collisions if event weighting is done only in 
| cos 9 1 the method provides the acceptance corrected 
asymmetry without applying any acceptance correc- 
tions for missing coverage in | cos 9\ . The | cos 9\ accep- 
tance fully cancels to first order. For proton-antiproton 
collisions, this is only true if the acceptance is the same 
for positive and negative values of | cos#. If it is not, 
then the small difference in the acceptance for posi- 
tive and negative | cos 9 (e.g. dead channels) needs to 
be accounted for in a MC acceptance correction. Ef- 
fects like a displaced beam vertex need to be corrected 
for using a MC simulation (for the proton-antiproton 
case). However, for proton-proton collisions, effects of 
the differences in the detector between forward and 
backward hemispheres or a displaced beam spot com- 
pletely cancel. This is because for proton-proton colli- 
sions the forward-backwards arc defined with respect 
to the rapidity of the dilepton, which can be either pos- 
itive or negative in the detector. A correction for the 
acceptance in \y\ still needs to applied in order to ex- 
tract the average integrated asymmetry over all values 
of \y\. This the method that should be used to extract 
the quark couplings from the average integrated asym- 
metry over all values of \y\. 

4. If the event weighting is be done in both |cos#| and 
| y |) then both the | cos#| acceptance and the \y \ accep- 
tance fully cancel. This is the method that should be 
used in a search for new Z' bosons. It should also be 
used for the extraction of the asymmetry in the pro- 
duction of top-antitop pairs, for which the acceptance 
is a complicated function of both | cos#| and \y\. 

5. Only small corrections need to be made are detector 
resolution and radiative smearing effects. Most of these 
correction are already included if an empirical fit Wi — 

f{y Unmeasured, Mu,measure.d) IS USed. 

6. Only very small corrections need to be made for QCD 
modification of the angular distribution. Most of these 
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corrections are already included if we use an empirical 
fit for h(8, Mu, Pt), as given in Appendix A. 

7. The method uses event weights which are functions 
of either \y\ or |cos#|. This is equivalent to weight- 
ing the data (in bins in \y\ or |cos#| ) by the error 
in the expected number of events. Therefore, this is 
equivalent to using the maximum likelihood method 
(which yields the smallest error). However, the imple- 
mentation of this technique is much easier and does 
not depend on a detailed modeling of the detector. 

8. If the event weighting is be done in both |cos0| and 
\y\) then both the |cos6>| acceptance and the \y\ ac- 
ceptance cancel. Therefore, for this case the weighted 
sums for Al, A2, Bl, B2, AB U and AB 2 from dif- 
ferent run conditions, or different final state leptons 
(muon, electrons) or different experiments (e.g. Dzcro 
and CDF or CMS and ATLAS) can be directly added 
to provide a combined result. This is important when 
there are only a few events in each detector (e.g. Drcll- 
Yan events at high mass or top-antitop events). 

8 Correction factors and systematic 
uncertainties in the weighting procedure 

For pp collisions, the acceptance for forward and backward 
events is equal because of symmetry. Therefore the func- 
tional dependence of the acceptance in | cos 8\ fully cancels 
(if the acceptance for positive and negative muons is the 
same). For pp collisions one needs to correct for a possible 
small difference in the detector performance between the 
proton and antiproton directions. If the event weighting 
is be done in both |cos#| and \y\ then both the |cos#| 
acceptance and the \y\ acceptance cancel. 

8.1 backgrounds and systematic errors 

The main experimental background is QCD dijet events. 
The QCD jet background is measured by statistically sep- 
arating isolated muons (or electrons) from muons (or elec- 
tron like objects) in jets on the basis of the transverse 
energy profile distributions in the calorimeter [T^] (e.g. 
isolation energy variables). 

In general QCD processes are mediated via spin 1 
gluon exchange and therefore have the same angular dis- 
tribution as Drell Yan events. If such is the case, the 
fractional QCD background is the same at all values of 
cos 8 and is the same for positive and negative values of 
cos 8. Therefore, in the expressions for the asymmetry, the 
QCD background cancels in the numerator. Corrections 
for this background only increases the level of the denom- 
inator. Therefore, a single multiplicative factor equal to 
1/(1 — /) (where f is the fraction of QCD background 
events) can be used to correct the extracted A /{, for QCD 
background. The same multiplicative correction factor for 
QCD background can be use for all three methods (simple 
count, fit to bins in | cos 8\ and event weighting technique). 
Possible deviations from this assumption can be included 
in the systematic error. 



Another background originates from electroweak (EW) 
processes (WW, WZ, M^+jets, t+t~, and ti). This back- 
ground is generally estimated from a Monte Carlo simu- 
lation. If the event weighting technique is used to extract 
the forward-backward asymmetry, then the event weight- 
ing technique can be used on Monte Carlo samples for 
the Drell- Yan signal and the electroweak background pro- 
cesses to determine the shift in from EW backgrounds. 
Note that the contribution from t + t~ events is very small 
(and at high mass Afi, for t + t~ is similar to Afi, for fifj, 
and e + e~ events). 

As is generally the case in particle physics experiments, 
the procedure needs to be tested on Monte Carlo simu- 
lated data to determine the size of any pulls from the pre- 
viously listed systematic effects. We note that this process 
to determine biases and pulls needs to be done for any 
procedure that is used to extract the forward-backward 
asymmetry for qq processes from data. Our procedure is 
designed to minimize these biases, but they still need to 
be determined from a full scale Monte Carlo simulation. 

In pp collisions at the LHC, there are also uncertain- 
ties from antiquark distribution functions that affect the 
misID probabilities (which are determined using a Monte 
Carlo simulations). The PDF uncertainties are smaller 
than the statistical errors for large dilepton mass. There- 
fore, this method is best applied to searches for new reso- 
nances or contact interactions at large mass. At large mass 
the leading order qq diagram dominates and the signal for 
new physics is local change in the asymmetry as a function 
of dilepton mass. 

However, for precision measurement of Standard Model 
parameters at the LHC (e.g. near the region of the Z bo- 
son) only weighting in | cos#| should be used. This is be- 
cause at lower dilepton mass, gluon-quark processes are a 
significant fraction of the cross section, and the dilution 
correction is much more sensitive to PDFs and higher or- 
der QCD terms. At low mass, it is better to compare the 
measured asymmetry (using |cos^| weighting) to a com- 
plete NLO or NNLO QCD prediction. 

9 Summary 

We have shown that a simple event weighting technique 
can be used to reduce the statistical error on the extracted 
qq forward-backward asymmetry from Drell- Yan events in 
pp and pp collisions. In addition to reducing the statistical 
error, the event weighting technique is much less sensitive 
to details of the experimental acceptance. 

The asymmetry data for a dilepton mass above 300 
GeV/c 2 (for both pp and pp collisions) can be used to 
search for new Z' bosons. 

10 Appendix A - Derivation of the 
perturbative formula for the angular 
coefficients. 

We now show that for the qq — > j* /Z — > i + £~ annihila- 
tion diagram the perturbative calculation indicates that 
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the angular coefficients originates from a kinematic effect 
which is a simple function of the angle of the proton and 
antiproton beams in the Collins-Soper frame. That angle 
is also a simple function of the transverse momentum Pt 
and mass Mm of the dilepton pair. 

We start with the case of the emission of a single hard 
gluon by one of the partons. In this case, the dilepton 
has a finite value of Pt but we cannot determine if the 
gluon was emitted by the parton in the proton, or by the 
parton in the antiproton. All we know is that the dilepton 
pair has a finite transverse momentum in the laboratory 
system which is equal to the transverse momentum of the 
gluon. 

If the gluon is emitted by the parton in the antiproton, 
then the parton level angular distribution (given below in 
cos 9') is with respect to the direction of the proton (which 
defines the quark-antiquark axis. If the gluon is emitted 
by the parton in the proton, then the parton level angu- 
lar distribution (given below in cos 9') is with respect to 
the direction of the antiproton (which defines the quark- 
antiquark axis) . The Collins-Soper frame splits this differ- 
ence and takes the average direction (which happens to 
be the same as the laboratory z axis for colliding beam 
experiments). 

In the derivation below, all vectors are defined in the 
Collins-Soper frame. We define unit vector k along the z 
axis, unit vector i along the x axis of the Collins-Soper 
frame (which is defined to be along the initial direction 
of the transverse momentum of the dilepton pair in the 
laboratory system). We also define a perpendicular unit 
vector j along the y axis. Here 8 is the angle with respect 
to the z axis, and <p is defined in the x — y plane. 

In the Collins-Soper frame, the Drell-Yan pair is sta- 
tionary and proton and antiproton move towards the dilep- 
ton pair with an angle £ with respect to the z axis. In the 
Collins-Soper, we see an antiproton moving from the left 
at an angle £ with respect to the z axis (in the x — z plane), 
and a proton moving from the right with an angle— £ with 
respect to the z axis. The proton and antiproton form an 
isosceles triangle with the two beam particles (which have 
equal momentum p) intersecting at an angle 2a — -k — 2£. 
Here cos a = sin £. It can be easily shown that 



sin 2 £ 



cos £ 



P 2 



P 2 T 



M/i 



(33) 



(P| 



Mj t yn 



In the Collins-Soper frame, we define a unit vector Z' 
along the direction of one of the beam particles (e.g. the 
proton). We define a unit vector r along the direction of 
Drell-Yan pair decay positron. We then obtain the follow- 
ing relations. 

Z' = cos£ k+ sm£ i (34) 
r = sin 9 cos <fi i + sin9 sincj) j + cos 9 k 
cos 6' = r Z' 

= cos £ cos 9 + sin £ sin 9 cos <j> 



For the case for which the gluon was emitted by the par- 
ton in the antiproton, the angular distribution of the decay 
positron with respect to the Z' proton axis (which defines 
the quark-antiquark axis) is given by the simple expres- 
sion. 



da 



= A 



(1 



?') + ^cos( 
A 



(35) 



e?cos(6>) 

(if the gluon was emitted by the proton, then the above 
distribution is with respect to the direction of the an- 
tiproton which defines the quark-antiquark axis). There- 
fore, the distributions for the proton direction and for 
the antiproton direction should be added. Each should be 
weighted by the probability that the quark emitted the 
gluon P q (xi), or the probability that the antiquark emit- 
ted a gluon, Pq(x2)- Here P q (xi)+Pq{x2)=1~ For proton- 
antiproton collisions, the two probabilities should be equal. 
For proton-proton collisions, they could be different. 

The proton direction is at an angle £. The antiproton 
direction angle is — £. If we plug in the expression for cos 9' 
from equation [34] into equation [35} we get 

da 



dco${6') 
term a 



term a + term* 



(36) 



cos 2 9 cos 2 £ 



AP q ( Xl )[l 
sin 2 9 sin 2 £ cos 2 <j) -f 
2cos£ cos #sm£ sin( 
B' 

sm£ sin 8 cos ( 



cos ( 



cos£ cos 9 



terrriq = CP q {x2)[l + cos 2 



sin 2 £ cos 2 . 



cos 2 £ 



2cos£ cos 9sin£ sin 9 cos< 
B 1 

— cos£ cos 9 
A 

B' 



A 



sot£ sind cos < 



Adding term q and term q and defining AP 
Pq{x2)) ~ we get 

da . rj 9 „ 9 . 

>s 2 (9 cos 2 £ 



dcos(9) 



A \1 



(37) 



sin* W sin 2 £ cos 2 <p + 
AP 2 cos £ cos 9 sm£ sin 9 cos < 
B 1 



A 



cos£ cos 9 



B' 

+ AP— sinf sink 
A 

using trigonometric identities 

sin 29 = 2 sink 

COS 2 (/)= -(1 + 

cos2<t>) 



COS t 



cost 



(38) 



(39) 
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we rewrite equation |37| as: 

= A[l + cos 2 



dcos(9) 



(40) 



+ ^(l-Scos 2 ^ 

+ AP cos £ sin £ sin 29 cos < 

sin 2 £ . 2 /> 2 , 
+ — ^ sin 2 6> cos 2 

+ AP— sin£ sin 6* cos 



IT 
T 



cos £ cos 9] 



The general expression for the total differential cross sec- 
tion for dilepton pair production (e.g. Drell-Yan, Z's or 
W's) in proton-antiproton (or proton-proton) collisions is 
given by: 



<r = A[(l + cos 2 9) 
+ 2^o(l - 3cos 2 ( 



Ai sin 29 cos < 



1 



A 2 sin 2 6* cos 2(f) + A 3 sin 6* cos < 



A4 cos + A$ sin 20 sin 2(j> 
A e sin 26* sin <f> + A 7 sin sin i 



(41) 



where A4 = -r . Comparing equation 40 to equation 
get: 



41 



we 



^0 
A 1 
A 2 
A 3 



sin 2 £= 152 



P 2 

T 



P't + M 2 U 



(42) 



AP cosf sin£ = AP _ , 



P T M t , 



= sin 2 £ = 



P 2 



Ml 



Pi 



B' B' 
AP— sin£ = AP— — ^5 ' „ N1 ,„ 



A 4 = T = 

P 
~B' 

A 5 



B _ B' B' M u 

A ~ T C ° S ^ ~ T(P 2 +M 2 )i/2 



(P| + M 2 ,) 1/2 
A 6 = A 7 = 



Where P is the cos 9 angular coefficient for quarks bound 
in the proton, and B' is the cos# angular coefficient for 
free quarks. As far as we know, this is the first closed form 
derivation of the QCD corrections to the forward back- 
ward asymmetry (i.e. the expression above for A4 and -Jv). 
The above perturbative formulae for Aq and A 2 for the 
annihilation diagram for Drell-Yan production was first 
derived by Collins in 19 79 [13]. The perturbative expres- 
sions for A , Ai and A 2 and AP (but not for A 3 and A4) 
were also derived in a paper by Daniel Boer and Werner 
Vogelsang in 2006) [13] . The above expressions for Aq and 



A 2 were shown to hold for the quark-antiquark annihila- 
tion diagram for all orders (re-summation) by Berger et. 
al. in 2007Q3]. 

In conclusion, if the Zj DY cross section is dominated 
by quark-antiquark annihilation (as expected) the pertur- 
bative calculation gives predictions for Aq (= A 2 ) and A4 
as the only surviving non zero angular coefficients. The 
simple QCD predictions for the angular coefficients are: 



h(9,M u ,P T ) = iA (l-3cos 2 
A 



(43) 



sin 2 £ 



A4 

B 
B 1 

cos £ 



sin 2 £ 
P 2 

T 



P 2 T 

B 

A = 



- M 2 

lvl u 

&_ 
~A 



cos£ 



cos£ 



M u 



(P 2 +M2)l/2 



The angular distributions predicted by standard Monte 
Carlo generators such as PHYTIA[9], ResBos re-summation [15 
and Vector Boson Production and Decay in Hadron Col- 
lisions Resummation VBP[16 for Z boson production at 
the Tevatron are well described by the above formulae. 

On average, both QCD correction terms are very small. 
For Drell-Yan events the average ratio Pt/Mu is about 
0.13 (e.g. 12 GeV P T at the Z mass of 90 GeV). This im- 
plies that on average < A (P T ) >= 0.018, and < -§t(Pt >= 
0.991. 



10.0.1 Including the effect of photon radiation by final state 
leptons 

Here, we derive an approximate formula for including the 
effect of final state photons. The more exact way is to run 
the Monte Carlo, with and without photon radiation and 
find a correction to the final result, empirical correction 
radiation. 

We assume that final state radiation smears the mea- 
sured angle of the final state electrons by an RMS angle 
±<5 in the 9 direction, and the same angle ±8 in <j>. The 9 
smearing is equivalent of taking the average of the above 
expressions with £ + 5 and £ — 5. The <j) smearing is equiv- 
alent of taking the average of the above expressions with 
2<f)+S and 2<fi—5. Therefore, the effect of final state photon 
radiation can be approximated by: 



A™ d = i[sin 2 (e + 5)+sin 2 (e-<5)] 

= sin 2 (£) cos 2 6 + cos 2 £! sin 2 S 
A r 2 ad = [sin 2 £ cos 2 6 + cos 2 £ sin 2 5} 

x [- [cos(2</> + 5) + cos(20 + 5)} 

= [sin 2 £ cos 2 5 + cos 2 £ sin 2 5] cos 6 



(44) 
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\ rad 



A r ad cos S 
[cos(C 



1 

"12 

= fracB'A cos £! cost) 



The above formula approximates the effects of photon 
radiation using one effective parameter S. That parame- 
ter depends on the specific experimental cuts and chan- 
nels. For example, if the final state is composed of two 
electrons, some of the photons are merged in with the 
electron shower in the calorimeter and help restore part 
of the original electron energy and direction. If the final 
state is a dimuon, the radiated photons are not included 
in the muon energy. Therefore S is only an effective pa- 
rameter that needs to be determined from fits to a Monte 
Carlo simulation. Probably, the simplest way is to include 
the corrections for the effect of photon radiation as part 
of an overall MC radiative correction. Alternatively, we 
could extract 6 from a fit to the MC and parametrize the 
QCD+QED corrections to the angular distribution by: 



1 



h(6,M u ,P T ) = - A 



rad 1 



rad 



A 



j^rad 
j^rad 

B> 
sin 2 £ 



cos £ 



1 - 3 cos 2 6) 
2 x , „„„2 . 



(45) 



sin £ cos S + cos £ sin S 
cos£ COS 6 



A A 
cos£ cos 8 

p2 



P 2 T 



-Mm 



(P 2 + M 2 ,) 1/2 



For 5 of order 0.1, the effect of QED photons is of the 
same magnitude as the effect of the QCD corrections. 
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